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INTRODUCTION 


In  millimeter  wave  devices,  a  solenoidal  magnetic  field  is  usually  employed  to  ensure  charged 
particle  beam  focussing.1' 4  The  magnetic  field  is  often  provided  by  electrical  solenoids  that  tend  to 
be  heavy  and  bulky  and  require  cumbersome  power  supplies.  The  reduction  of  size  and  weight  of 
millimeter  wave  devices  has  been  increasingly  emphasized,  especially  for  military  applications. 
Permanent  magnet  solenoids  that  employ  radially  oriented  cladding  magnets  to  confine  flux 
produced  by  axially  oriented  supply  magnets  have  been  the  subject  of  recent  studies.5'8  In  addition, 
such  clad  magnetic  structures  have  found  applications  in  many  other  devices;  for  example,  magnetic 
resonance  imagers,9  an  electron  filter  magnet10  and  an  X-ray/UV  image  sensor.11 

The  nature  of  the  application  dictates  the  strength  of  the  desired  field  and  the  dimensions  of  the 
work  space.  These,  in  turn,  determine  the  dimensions  of  the  solenoid  and  the  configurations  of  its 
components.  Recent  applications  have  required  the  ability  to  design  many  solenoids  of  varying 
specifications,  through  the  use  of  a  standard  set  of  basic  components. 

Accordingly,  it  was  decided  to  compile  a  catalog  of  permanent  magnet  solenoids  to  be  used  as  a 
handy  reference  by  the  user  in  need  of  such  a  device.  It  provides  basic  information  with  regards  to 
mass,  volume  and  workspace  fields  and  affords  the  user  a  basis  for  assessing,  a  priori,  the  pros  and 
cons  of  the  permanent  magnet  design  over  its  electrical  counterpart. 

THEORY 

Figure  1  shows  a  clad  cylindrical  permanent  magnet  structure,  where  the  flux  is  provided  by  an 
axially  magnetized  permanent  magnet  shell,  a,  encompassing  the  working  space,  W .  The  shell  is 


capped  at  its  ends  by  iron  discs,  b,  which  guide  all  of  the  magnet’s  flux  Om  into  W.  Figure  2 
summarizes  the  calculation  of  the  supply  magnet  thickness  from  the  desired  field  and  the  dimensions 
of  the  working  space.  From  Ampere’s  Law  and  the  assumption  of  perfect  flux  confinement,  we  have 
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Figure  2:  Configuration  of  Supply  Magnet  and  Pole  Pieces  in  the  Workspace 


®m=®w 

H m  =  Hw  (!) 

where  Oyy  is  the  flux  in  the  working  space,  Hm  is  the  field  within  the  supply  magnet  and  Hw  is 
the  field  in  the  working  space.  Equating  the  flux  inside  the  supply  magnet  to  the  flux  within  the 
working  space,  we  have 


Byn^m  ~  H-w^w 


(2) 


where 


Byn  ~  VrH m  By  (3) 

and  Am  and  Aw  are  the  cross-sectional  areas  of  the  supply  magnet  and  of  the  working  space, 
respectively.  fJLr  denotes  the  slope  of  the  demagnetization  curve  of  the  magnet  material  used  and  is 
approximately  1  for  rigid  permanent  magnets.  Substituting  for  Bm  in  Eq.  (2)  and  solving  for  Am 
we  obtain 


Am  ~  Aw/(pr  Br/\HW\)  (4) 

Equation  (4)  is  written  in  the  form  shown  to  obviate  the  need  to  include  the  sign  of  Hm,  since 
Hm  <  0.  The  outer  radius  of  the  supply  magnet  is  determined  by  Am  through  the  expression 


where  rQ  and  r(  are  the  outer  and  inner  radii  of  the  supply  magnet,  respectively. 
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Figure  3:  Cladding  Magnet  Construction 


To  keep  the  flux  confined  to  Ayy  we  need  a  tapered,  radially  oriented  cladding  magnet  encasing 
the  annular  longitudinal  supply  magnet.  The  supply  magnet  flux  must  be  confined  by  a  cladding 
magnet  of  unknown  thickness  d(x),  shown  in  Fig.  3.  The  thickness  of  the  cladding  magnet  is 
determined  by  the  circuital  form  of  Ampere’s  Law,  which  states  that,  in  the  absence  of  linking 
currents,  jH •  dl  =  0  around  a  closed  loop,  where  dl  is  an  infinitesimal  length  in  the  direction  of 
the  closed  path  and  H  is  the  magnetic  field  along  the  same  path.  Therefore, 

0  =  f  Hdl  =  Hwx  +  Hcd  +  Hext(BO)  (6) 

where  Hq  is  the  radial  field  in  the  cladding  magnet  and  Hext  the  field  outside,  which  is  zero  by 
hypothesis.  Hence, 

d  —  —  Hmx/Hc  (7) 

where  Hc  is  just  the  magnet  coercivity  and  d  is  the  cladding  thickness  at  any  point  x.  Equation 

(7)  yields  the  cladding  thickness  at  any  point  along  the  device  for  0  <  X  <  L  where  L  is  the  total 
length  of  the  structure. 

Next,  to  prevent  end  losses  we  must  have  pole  pieces.  We  find  the  thickness  of  the  pole  pieces 
by  equating  the  flux  in  the  working  space  to  the  flux,  d?p,  traversing  the  circular  band  that  forms 
the  circumference  of  the  pole  piece.  The  result  is  : 

Q  P  =  HyyAyy  =  'IjQyyt pB p 
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(8) 


tp  —  HyjAyj  I'lTtTyyBp 

where  rw  is  the  radius  of  the  working  space  and  tp  is  the  thickness  of  the  pole  piece.  If  the  pole 
piece  is  to  be  unsaturated.  Bp  must  be  no  larger  than  Bs,  the  saturation  induction  of  the  pole  piece 

material.  In  our  case,  Bs  is  the  saturation  induction  of  iron  (20  kG).  To  build  in  a  safety  factor  to 
prevent  local  saturation,  we  multiply  the  resulting  minimal  thickness  by  2  and  obtain 

tp  =  HwAw/xrwBs  (9) 

With  the  determination  of  these  dimensions,  the  design  of  the  solenoid  can  be  completed. 

To  get  maximum  confinement,  we  must  also  include  in  the  design  disc-shaped  bucking  magnets 
and  ring-shaped  comer  magnets.  The  thickness  of  the  end  cladding  is  the  same  as  the  maximum 
cladding  thickness,  tc  (max) ,  calculated  for  x  =  L  in  Eq.  (7) 

FINITE  ELEMENT  ANALYSIS 

To  cover  a  wide  range  of  requirements,  solenoids  were  modeled  for  working  space  fields  ranging 
from  500G  to  6500G,  and  aspect  ratios  ranging  from  0.5  to  10.0.  Here,  aspect  ratio  is  defined  as  the 
ratio  of  length  to  diameter  of  the  working  space.  All  the  solenoids  modeled  had  a  working  space 
diameter  of  2.0  cm  and  were  made  of  Nd-Fe-B.  A  remanence  of  12  kG  for  all  permanent  magnet 
materials  and  density  of  8.3  g/cc  for  all  materials  was  assumed  in  the  entire  design. 

Finite  element  analysis  was  performed  using  2D  magnetic  analysis  software,  version  13.1,  from 
the  MacNeal  Schwendler  Corporation.  A  typical  plot  of  the  magnetic  flux  is  shown  in  Fig.  4. 


- I .  t _ E_: _ l|. 

Figure  4:  Magnetic  Flux  Distribution  in  a  Permanent  Magnet  Solenoid 
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RESULTS 


The  flux  plot  of  Fig.  4  shows  excellent  flux  confinement  by  the  radially  oriented  cladding 
magnets  and  the  disc-shaped  bucking  magnets.  There  is  some  leakage  of  flux  at  the  comers  near  the 
ring-shaped  bucking  magnets.  This  issue  has  been  dealt  with  elsewhere12  and  can  be  reduced  with  a 
modification  to  the  basic  solenoid  if  the  leakage  presents  a  problem  to  instrumentation  placed  in  that 
region. 

Figures  5-11  show  the  field  profiles  along  the  axis  of  the  solenoids,  sequentially,  for  fields 
ranging  from  500  gauss  to  6500  gauss  in  increments  of  1000  gauss.  The  field  estimated  by  the  finite 
element  calculation  is  indeed,  on  the  average,  the  field  expected  from  the  design.  It  would  seem 
from  a  quick  glance  at  the  field  profiles,  however,  that  there  is  a  significant  deviation  from  the 
desired  field  along  the  axis.  The  fields  are  somewhat  higher  near  the  center  of  the  solenoid  and  drop 
off  considerably  towards  the  ends  of  the  axis.  For  small  aspect  ratios,  i.e.  <  1,  the  percentage  of 
error  ranges  anywhere  from  1.5%  for  small  fields  to  3.0%  for  the  larger  fields.  For  larger  aspect 
ratios  these  percentages  of  error  are  less  impressive,  starting  from  a  mere  5.2%  for  a  500G  solenoid 
and  rising  to  to  38.5%  for  our  longest  design,  an  aspect  ratio  of  10,  for  the  highest  field  in  this  study, 
6500  G.  If  such  deviations  in  polar  axis  fields  are  deemed  too  large,  remedies  are  available  to  make 
the  field  profile  more  uniform.12 


Figure  5:  Axial  Fields  for  a  500  gauss  Solenoid 
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Axial  field  in  gauss  Axial  field  in  gauss 


4500 


Distance  along  the  axis  in  centimeters 


Figure  8:  Axial  Fields  for  a  3500  gauss  Solenoid 


Figure  9:  Axial  Fields  for  a  4500  gauss  Solenoid 


Axial  field  in  gauss  Axial  field  in  gauss 


Figure  12:  Percentage  of  Error  in  the  Obtained  Field  as  a  Function  of  Aspect  Ratio 


Figure  13:  Solenoid  Mass  as  a  Function  of  Desired  Field 


Figure  12  summarizes  the  difference  in  the  required  fields  and  the  peak  fields  at  the  center  of  the 
solenoid,  expressed  as  a  percentage  of  error,  as  a  function  of  aspect  ratio.  If  the  length  of  the 
working  space  is  small  in  comparison  with  the  inner  radius  of  the  solenoid,  annular  confinement  of 
flux  is  favored. 

CONCLUSION 

A  permanent  magnet  solenoid  can  be  designed  accurately  for  any  desired  field.  If  field 
uniformity  is  the  criterion  for  evaluating  the  success  of  the  design,  the  permanent  magnet  solenoid  is 
best  for  small  fields  and  aspect  ratios.  As  the  requirements  shift  towards  higher  working  space 
fields,  one  must  sacrifice  uniformity  in  field  profile,  particularly  if  the  demand  is  for  a  working  space 
that  is  several  fold  longer  than  it  is  wide.  Figure  13  is  the  designer’s  guide  to  the  increase  in  the 
mass  of  the  solenoid  as  one  demands  a  larger  field. 

It  would  appear  from  studying  the  magnetic  field  profiles  in  Figs.  5-11  as  well  as  the  summary 
presented  in  Fig.  12,  that  field  uniformity  in  the  interior  of  the  solenoid  can  be  rather  poor, 
particularly  for  high  aspect  ratios.  This  arises  due  to  imperfect  cladding  at  the  ends.  Many 
applications  find  this  field  gradient  acceptable  while  others  are  stringent  in  their  demands  for  a 
uniform  field.  For  the  latter,  the  permanent  magnet  solenoid  is  still  a  viable  alternative.  A  simple 
modification,  discussed  in  detail  elsewhere,9’12  can  improve  the  field  profile  dramatically.  This  is 


Figure  14:  Axial  magnetic  field  profile  in  a  cladding  compensated  solenoid 
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accomplished  by  altering  the  thickness  of  the  cladding  magnet.  It  can  be  done  geometrically.  The 
correction  consists  of  increasing  the  cladding  thickness  At  according  to  the  formula: 

A  t  =  (Hd-Ha)t/Ha  (10) 

where  Hrf  is  the  desired  field  and  Ha  is  the  actual  field.  Figure  14  shows  the  axial  field  as  a 
function  of  distance  from  one  end  in  a  sample  solenoid.  The  much  improved  field  uniformity 
obtained  after  a  simple  correction  to  the  cladding  slope  is  made  is  also  included.  The  deviation  from 
the  average  field  in  the  unaltered  solenoid  is  3.0%  and  is  reduced  to  less  than  0.01%  upon  geometric 
modification  of  the  cladding. 

It  is  also  possible  to  vary  the  cladding  parametrically.3  Br,  the  remanence  of  the  cladding 
magnets,  can  be  gradually  increased  in  the  manner  shown  below: 

ABr  =  HaAt/t  (11) 

where  At  is  the  increment  in  cladding  thickness  t  determined  by  Eq.  (10). 
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